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Abstract. Laurent polynomials related to the Hahn-Exton g-Bessel function, which are q- 
analogues of the Lommel polynomials, have been introduced by Koelink and Swarttouw. The 
explicit strong moment functional with respect to which the Laurent gr-Lommel polynomials 
are orthogonal is given. The strong moment functional gives rise to two positive definite 
moment functionals. For the corresponding sets of orthogonal polynomials the orthogonal- 
ity measure is determined using the three-term recurrence relation as a starting point. The 
relation between Chebyshev polynomials of the second kind and the Laurent q-Lommel poly- 
nomials and related functions is used to obtain estimates for the latter. 

> 

£^ ■ 1. Introduction and motivation 

o : 

l/-) . The Lommel polynomials are orthogonal polynomials closely related to the Bessel func- 

tion. Although the Lommel polynomials have a representation involving a hypergeometric 
2i 7 3-series, they do not fit into Askey's scheme of hypergeometric orthogonal polynomi- 
als. The reason for this is that the orthogonality measure for the Lommel polynomials is 
supported on the set consisting of one over the zeros of a Bessel function, which are not 
>• . explicitly known in general. So there is no Rodrigues formula or difference equation for 
the Lommel polynomials. 

The Bessel function J v {z) of order v and argument z is given by the absolutely conver- 
gent series expansion 

k\r(v + k + i) 
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The properties of this special function are well understood, see e.g. the book on Bessel 
functions by Watson [22]. A simple recurrence relation for the Bessel functions is, cf. [22, 
§3-2(1)], 

(1.2) J u+1 (z) = —J u (z)-Ju-i(z). 

z 

From iteration of (1.2) we see that we can express J u+m (z) in terms of J v {z) and J v _\{z) 
and the coefficients of J v [z) and J v _\{z) are polynomials in z~ x . This was first observed 
by Lommel in 1871. Explicitly, we have, cf. Watson [22, §9.6], 

(1-3) J u+m (z) = h mu (-)J v (z) - h m _ l l/+1 (-)J l/ _ 1 (z), 

z z 

where h m ^ v (z) are the Lommel polynomials, which are also known as associated Lommel 
polynomials. The Lommel polynomials satisfy the three-term recurrence relation 

(1.4) h m +i iV (z) = 2z(m + v)h m , v (z)-h m -i jV (z), h- 1>v (z) = 0, h , u {z) = 1. 

Favard's theorem, cf. Chihara [7, Ch. II, thm. 6.4], implies that the Lommel polynomials 
are orthogonal polynomials with respect to a positive weight function for v > 0. The 
explicit orthogonality relations are, cf. Chihara [7, Ch. VI, §6], Dickinson [9], Dickinson, 
Pollak and Wannier [10], Ismail [15], Schwartz [19], 



(L5) £ or 1 ) 2 hm,v Or 1 ) v " Or 1 ) " 2 (-+-)' 

where v > —1, are the positive zeros of the Bessel function J v (z) numbered increasingly, 
cf. Watson [22, Ch. 15]. 

Another relation between the Lommel polynomials and the Bessel function is given by 
Hurwitz's asymptotic formula, cf. [22, 9.65(1)], 

1.6 — ' >J^_i (-), m^oo. 

T(u + m) \zJ 

For the Bessel function (1.1) there exist several (/-analogues. The oldest (/-analogues for 
the Bessel function were introduced by Jackson in a series of papers in 1903-1905, see the 
references in [15]. For the Jackson (/-Bessel function Ismail [15] introduced the associated 
g-Lommel polynomials, which turned out to satisfy an orthogonality relation similar to 
(1.5), but now involving the zeros of the Jackson (/-Bessel function. Ismail used these 
(/-Lommel polynomials to prove that the zeros of the Jackson (/-Bessel functions behave 
similarly as the zeros of the Bessel function. 

A more recent (/-analogue of the Bessel function has been introduced by Hahn in a special 
case and by Exton in full generality, see the references in Koornwinder and Swarttouw [18]. 
The zeros of the Hahn-Exton (/-Bessel function and several associated (/-analogues of the 
Lommel polynomial have been studied by Koelink and Swarttouw [17]. The zeros of 
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the Hahn-Exton (/-Bessel function behave in a similar fashion as the zeros of the Bessel 
function. In that paper [17] a (/-analogue of the Lommel polynomials was introduced. 
However, this (/-analogue of the Lommel polynomial is no longer a polynomial, but a 
Laurent polynomial. One of the goals of this paper is to give an explicit orthogonality 
measure for these orthogonal Laurent (/-Lommel polynomials. 

The Laurent (/-Lommel polynomials are defined by, cf. [17, prop. 4.3 with R rrijl/ (z~ 1 ; q) = 
h m ,v(z;q)], 

(1.7) h m+ltV (x; <?) = (- + x(l - q l ' +m ))h m ^{x] q) - h m _ x ^{x\ q), 

with initial conditions h-i yU (x;q) = 0, ho :ll (x;q) = 1. A second independent solution of 
(1.7) is given by /i m _i )I/ _ ) _i(x; q). Note that taking the limit q f 1 in (1.7) after replacing 
x by 2z/(l — q) gives (1.4). The Laurent (/-Lommel polynomials originate from a relation 
similar to (1.3), see proposition 3.1. 

The explicit orthogonality relations for the Laurent (/-Lommel polynomials h m ^{x\q) 
defined in (1.7) is derived in §3. The method of proof is based on the existence of asymp- 
totically well-behaved solutions of (1.7) reminiscent of J u+m (x), cf. (1.3). The method 
of Dickinson [9] to prove (1.5) can then be adapted to our situation. The orthogonality 
measure gives rise to a strong moment functional £, i.e. a functional on the space of Lau- 
rent polynomials so that all moments £(x n ), n G Z, exist. From C we obtain two moment 
functionals £±, as considered in e.g., [7, Ch. 1], by putting C + (x n ) = C(x n ), n G Z + , and 
C-(x n ) = —C(x~ 2 ~ n ), n G Z_|_. (The 2 has to do with the fact that all moment func- 
tionals are symmetric.) It turns out that both C + and £_ are positive definite moment 
functionals. 

The orthogonal polynomials for C + are (/-analogues of the Lommel polynomials and 
the support of orthogonality measure consists of zero and one over the zeros of a Hahn- 
Exton (/-Bessel function. This is worked out in detail in §4, where we use Dickinson's 
method [9] once more. In §5 we study the orthogonal polynomials for We give explicit 
expressions for these polynomials in terms of Al-Salam-Chihara polynomials, which can be 
used to determine the asymptotic behaviour as the degree tends to infinity. The asymptotic 
behaviour is expressed in terms of a function j v (x;q) closely related to the Hahn-Exton 
(/-Bessel function. Since we can do this for the associated polynomials as well, we have 
the Stieltjes transform of the orthogonality measure from which the orthogonality follows. 
Using the results of §5 we can simplify the expression for the strong moment functional C 
using a Wronskian type formula. This is done in §6. 

For q = 0, or for v — > oo, we see that U m ((x+x _1 )/2) , where £7 m denotes the Chebyshev 
polynomial of the second kind, satisfies (1.7) with the same initial conditions. So we can 
view the Laurent (/-Lommel polynomials h mjU (x;q) as a perturbation of the Chebyshev 
polynomials. This point of view allows us to obtain estimates for the Laurent (/-Lommel 
polynomials, the Hahn-Exton (/-Bessel function and the related function j l/ (x;q). This is 
done in §7. 

Finally, in §2 we show that the general theory of orthogonal Laurent polynomials 
presents us with an existence theorem for the strong moment functional C. We also state 
a result concerning the zeros of the Laurent (/-Lommel polynomials. 
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To end this introduction we briefly recall the notation for basic (or (/)-hypergeometric 
series. We follow the standard notation of Gasper and Rahman [11, Ch. 1]. We take 
< q < 1 for the rest of the paper. A (/-shifted factorial is a product defined by 



k-1 



(a;q) k =Y[(l-aq l ), a G C, k G Z+, 



i=0 



where the empty product equals 1 by definition. Since < q < 1 we can take k — > oo to 
get liuik^ QO (a;q)k = (a; <z)oo- A basic (or (/)-hypergeometric series is 



OO 



^ (V<l)k(h;q)k---(bs;q)k y ' 

For generic values of the parameters the radius of convergence of the series in (1.8) is 0, 1, 
oo corresponding to r > s + 1, r = s + 1, r < s + 1. 

2. Orthogonal Laurent polynomials 

In this section we apply some of the theory of orthogonal Laurent polynomials to the 
Laurent polynomials h mjl/ (x;q) to obtain the existence of a strong moment functional £, 
i.e. a linear functional on the space of Laurent polynomials for which the moments C(x m ) 
exist for all m G Z, for which the Laurent g-Lommel polynomials are orthogonal. We 
use the paper by Hendriksen and van Rossum [14] as the main reference for this section. 
The recurrence relation as in (2.1) has been generalised to a wider class of recurrence 
relations by Ismail and Masson [16] by replacing x in front of the V m -i jU (x) by (x — a m ), 
for which they prove a Favard type theorem. Specialisation to the case considered here 
yields the Favard type theorem contained in Hendriksen and van Rossum [14] . For further 
information concerning this section the reader may consult the introductory paper by 
Cochran and Cooper [8]. 

From the recurrence relation (1.7) it follows that h m ^{x\ q) is an even function for even 
m and an odd function for odd m. Consequently, x m h mjU (x;q) is a polynomial in x 2 , 
which we denote by V mjl ,(x 2 ) = x rn h rritU (x; q). For V m we obtain from (1.7) the recurrence 
relation 

(2.1) V m+1 , v (x) = (1 + x(l - q v+m ))VmAx) - xV m - ltV (x) 

with initial conditions V-\^(x) = 0, Vq jU (x) = 1, cf. [14, (2.2)]. The Favard-type theorem, 
cf. [14, thm. 1.1], implies that for the Laurent polynomials Q n (x), defined by 

Q 2 n(x) = x~ n V 2n , v (x) = h 2n ,v(Vx; q) 
Q2n+i(x) = x~ n ~ 1 V 2n +iA x ) = x~^h 2n+ljV (^/x;q) 
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there exists a strong moment functional C\ such that Ci{Q n Qm) = for n 7^ m. 

If we form the lacunary Laurent polynomials, cf. [14, (1.16)], we get the Laurent polyno- 
mials P2m(x) = h mtl y(x; g), P2 m +i(x) = x~ 1 h m ^(x; q). The lacunary Laurent polynomials 
are orthogonal with respect to the strong moment functional C defined by C(x 2n ) = £i(x n ), 
C(x 2n+1 ) = for n G Z, cf. [14, prop. III]. So the orthogonality relations for the even 
lacunary Laurent polynomials gives 

f == 0, 72 7^ 772, 

(2.2) C(h niV {x;q)h m , v {x;q))\ ' 

[ 7^ 0, n = m. 

But we also have the orthogonality for the odd lacunary Laurent polynomials, 

(2.3) C(x~ 2 h njl/ (x;q)h m:l/ (x;q))l °' n ^ m " 

The space A n of Laurent polynomials of the form Ylp=- n c p xP * s + 1) -dimensional, 
n G Z_|_. The Laurent polynomials h m ^(x] q), m = 0, 1, . . . , n, form a (n + 1) -dimensional 
subspace of A n . Moreover, they form an orthogonal basis for this subspace with re- 
spect to C. Equation (2.3) states that this orthogonal basis can be complemented with 
x~ 1 hm,v(x', q), m> = 0, 1, . . . , n — 1, to give an orthogonal basis of A n with respect to C. 
C(x~ 1 h mjU (x; q)h ny v(x; q)) = is immediate from C(x 2p+1 ) = 0. 

Remark 2.1. For orthonormal polynomials the three-term recurrence relation can be used 
to prove that the zeros of the orthonormal polynomials correspond precisely to the eigenval- 
ues of a truncated Jacobi matrix. A similar approach can be used here. Define coefficients 
by 

n+l 

(2.4) xV n>v (x) = ^2c n}k V k ,v(x), 

fc=0 

then the matrix H n = (cij)o<i,j< n -i is Hessenberg matrix, i.e. Cij = for i + 1 < j. 
Using (2.4) in (2.1) gives recurrence relations for the matrix elements c^j, which can be 
solved to give 

i 

if k = n + 1, 
if < k < n, 
if k = 0. 



(2.5) c n , k = < 



{q v ;q)n+i 
-l 



(q u ;q)n+i 

Note that each row sum of H n , except the last, equals zero. 

Introduce the vector w n (x) = (Vo jU (x), V\ jV (x), . . . , V n -i^{xj) , then we see from (2.1) 
that H n w n (x) = xw n (x) if V Hjl/ (x) =0. So a zero x of V n>u implies that H n has an eigen- 
vector for the eigenvalue x. It is also possible to prove that an eigenvalue x of H n implies 
that V n ^ v (x) = 0, which can be proved by showing that the characteristic polynomial of 
H n times the normalisation constant (— l) n (q u ; q) n satisfies (2.1). So we conclude that the 
zeros of V n ^ v (x), and hence the zeros of the Laurent g-Lommel polynomials h n ^(x; q), are 
completely determined by the spectrum of the Hessenberg matrix H n . 
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3. Minimal solutions and orthogonality relations 

In this section we give an explicit formulation for the strong moment functional L 
introduced in the previous section. We describe C in terms of contour integrals where the 
integrands depend on the Hahn-Exton y-Bessel function and on a function closely related 
to the Hahn-Exton g-Bessel function. These functions give rise to two other solutions of 
the recurrence relation (1.7), but now with prescribed behaviour for m — > oo. The proof of 
orthogonality of the Laurent g-Lommel polynomials for C uses a method already introduced 
by Dickinson [9] to prove the orthogonality relations (1.5) for the Lommel polynomials. 

Using a generating function argument the following explicit expressions for the Laurent 
g-Lommel polynomials have been derived in [17, (4.23)] from the recurrence relation (1.7) 



n+l 



(3.1) 

m / „n—m n n+l 

(6.2) — 2_^x 2 <fi I ,q,q 

n=0 ^ ^ 

The Hahn-Exton y-Bessel function is defined by 

(3.3) J u (x; q) = x i<pi I +1 ; q, qx 
and then the following (/-analogue of Hurwitz's formula (1.6) holds 

(3.4) lim x- m h m ^{x-q) = x v ~ x J,_! (-; q) , |x| > 1. 

Relation (3.4) has been proved formally in Koelink and Swarttouw [17, (4.24)] from (3.1), 
but it follows from their proof that it is valid only for |x| > 1. 

In order to state the asymptotic behaviour of the Laurent y-Lommel polynomials inside 
the circle we introduce the function 
(3.5) 

j v (x\ q) = x u {qx 2 ] g)oo i^i (0; qx 2 ] q, q v+l x 2 ) = x u (q u+1 x 2 ; q)^ i^i(0; q v+l x 2 ] q, qx 2 ), 

where we use (x; q)^ i<^i(0; x; q, y) = (y; q)^ i<^i(0; y; q, x), cf. [18, (2.3)]. This function is 
related to the Hahn-Exton y-Bessel function in the following way 



fi=u+2 In x/ In q 



x "ju(x;q) = (q;q) 00 (x M J M (x;g) 
Now we can use (3.2) to obtain 

x m h m ^x; q) = J2 x2m ~ 2U ( q ' * ; 9, q" +m ~ n 

■n=n V 9 / 



n=0 

™L / -n m-n+1 

= }_^x 2 <pi ( ;9,9 

n=0 ^ ^ 
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and by dominated convergence we obtain 



00 / ~ n n \ 

lim x m h m , v {x- q) = T] x 2n 2V ?i P ' ; q, q u+n 



00 



(9; 9)1(9; 9)1, _ 

where the last equality follows from interchanging the summations, which is allowed for 
\x\ < 1. The inner sum can be written as 

J2(Q~ P ~ l -i q)ix 2{p+l) q l(p+l) = x 2l (-l) l q^ l -V q' 1 )^ 

p=0 p=0 

= ^(- 1 )'9*'<«'(9;9) 1 E^^- 2 '(-l)'9i'<-)^, 

by the (/-binomial theorem, cf. [11, (1.3.2)]. This leads to the result 

1 / \ x x ~ v 

(3.6) lim x m h m ^{x]q) = - ^ 1^1 ( 2 ; 9, <?^ 2 = 7-3 — r— |x| < 1. 

m^oo 1 — x z \qx^ J (x ,q)oo 

Proposition 3.1. The functions J u + m (x~ ;q) and j u + m (x; q) satisfy the recurrence rela- 
tion (1.7). Moreover, 

J v + m {x~ x ; q) = h m , u (x; q)J u {x~ 1 \ q) - h m ^ u+1 (x\ q)J v _ 1 (x~ 1 \q), 
j u +m(x; q) = h m , u (x; q)j„{x; q) - h m - ljU +i(x; q)j v -i(x; q). 

Proof. Since h m ^{x;q) and h m -i^ + i(x; q) are linearly independent solutions of the re- 
currence relation (1.7), the last statement of the proposition implies the first. Also, if 
J l/ + m (x~ , q) and j l/ + rn (x;q) satisfy (1.7), then they must be a linear combination of 
hm,v{x', q) and h m -i tU +i(x; q) from which the second statement follows by considering the 
cases m = and m = —1. 

The last statement for J v + m has already been proved in [17, (4.12)], so it remains to 
consider j u + m . The second order (/-difference equation for the i</?i-series, or by taking a 
suitable limit in [17, (4.14)] in combination with (3.6), reveals that 

1 

jv+i(x;q) = (- + x(l - q v ))j v (x]q) -j u - 1 (x;q). 
Replacing v by v + m proves the statement. □ 

Remark, (i) The solutions J u+m (x~ 1 ; q) and j v+m (x; q) of (1.7) have the following asymp- 
totic behaviour for m — > 00 valid for x G C; 

,-2. 



lim x m+v J v+m (x- 1 ;q)= (g f \ q) 
(3.7) m ^°° (9; 9)oo 

lim x~ rn ~ u j u+m (x]q) = (qx 2 - 1 q) OCl . 



00 

00 
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Note that x ±m are solutions of (1.7) for m — > oo (or for q = 0). So the solutions 
J l/+m (x~ 1 ; q) and j u+m (x; q) behave as x^ m up to a factor independent of m as m — > oo. 
(ii) The functions J v s rrn {x~ 1 ] q) and < 7 l/ + m (x; g) are related to a minimal solution X m (x) of 
(2.1), i.e. X m (x) is a solution such that lim,™^^ X m (x) /V mjV (x) = 0, where Vm tU (x) is the 
polynomial solution of (2.1). Using the limit transitions (3.4) and (3.6) and the relations 
in proposition 3.1 we obtain 

' 3vWx;q)V mtV {x) - x^j v _ 1 {y/x;q)V m _ 1 ^ +1 {x) 

v , s = x* m j v+m (y/x;q), \x\<l, 

X m {x) — < i 

J v (l/y/x;q)V m ,v(x) - x»J v -i(l/y/x;q)V m -i,u+i(x) 

= x^ m J v+rn {l/^x;q), \x\>l. 



With the functions J u (x;q) and ju(x;q) and their relation with the Laurent g-Lommel 
polynomials described in proposition 3.1 at hand, we can give an explicit expression for 
the strong moment functional C. The proof we give is an adaption to the Laurent case of 
Dickinson's proof of the orthogonality (1.5) for the Lommel polynomials [9]. 

First we investigate the quotient of two Hahn-Exton g-Bessel functions. 

Lemma 3.2. For v > the following expansion holds around for n G Z + 

J v -\-ni,X i q) x "'" 



2k 

, I " I I 

J„-i(x;q) (q v ;q) n+1 j-f 



k=o 

with the coefficients c/- recursively defined by cq = 1 and 

(-l)M fc(fc+1) ^l)k- Pq \(k-p){k-p+l) 
(3 ' 8) ° k= (q» +n+1 ;q)k(q;q)k (q»;q)k(q;q)k- P 



Proof. From (3.3) we immediately get 

Ju+n(x;q) = x n+1 1 ^ 1 (0;q"+ n + 1 ;q,qx 2 ) 
J v -i(x;q) (q v ;q) n +i i¥>i(0; q v ; q, qx 2 ) 

so we have to solve for the coefficients Ck by comparing powers of x on both sides of 

- 2 ~ (_1)P Q |P(P+1) X 2 P ^ - (_ 1 )m g lm (m+l) a .2m 
2_^C k X („v.„\_(„.„\_ ~ 2^ (qv+n+l., 



k=0 



~ (q";q)p(q;q) P 



m=0 



L ;q)m(q;q)r 



from which the recurrence relation (3.8) for the coefficients c k is obtained. 
A rude estimate gives 



(-l) k qh k ( k + 1 ) 



<A = 



1 



(q";q)oo(q;q)c 
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for v > 0. The same estimate applies to the factor in front of c n on the right hand side of 

(3.8) , so that we obtain 

fc-i 

\c k \ <A + J2 A K\- 

n=0 

A discrete version of Gronwall's inequality, cf. e.g. [21, p. 440] 

k-l k-1 

(3.9) a k < A + ^2 d n a n , and A, a n , d n > =^ a n < A exp rf n ) , 

n=0 n=0 

yields |cfc| < Ae kA so that the series on the right hand side of the statement of the lemma 
is absolutely convergent for \x\ < e~ A / 2 . □ 

Choose < R < , where f\~ x denotes the smallest positive zero of J v -i{x;q), 
v > 0, cf. [17, §3]. Using lemma 3.2 we obtain for v > 0, m G Z and n G Z + 

(310) J_/ z m Jy+n{z- x ;g) dz _ / 0, m - n odd or m < n, 

2iriJ ]zl= i_ Ju-^z-^q) I m = n. 

Note that the coefficients c& of lemma 3.2 for n = are in fact the moments of the linear 
functional £+ defined by 

1 f JJz-^-iq) f 0, m 6 Z + odd, 

(3.11) £+(a: m ) = — i z m is \ ,1) dz= \ + 

2m J\ z \=i_ Jv-i{z l ;q) I Cm/2, m G Z + even. 

We return to this moment functional in section 4 and we calculate the corresponding 
orthogonal polynomials, which turn out to be g-analogues of the Lommel polynomials. 

The following lemma is the analogue of lemma 3.2 for the functions j u (x; q) instead of 
the Hahn-Exton g-Bessel function. 

Lemma 3.3. For v G R the following expansion holds around for n G Z + 

jv+n(x] q) _ x n+l J. J2,k 



with the coefficients dk recursively defined by do = 1 and 

k-i 

(3.12) d k = 2?i{q-\ 0; q; q, <f +"+*+*) - £ dp 2 ^i{q p -\ 0; q; q, ? "+ fc -*) 

Proof. The proof is completely analogous to the proof of lemma 3.2 and we only give the 
differences. Here use the expansion 

m (0; ,x 2 ; g, <f + V) = £ / *7 2 . ^ +1 >V fc 
(l-^ 2 ) (g;g) fc (a; 2 ;g) fc+ i 



fc=0 



= l~-U g 2 q (v+i)k x 2k (.q ,qn x 



oc 



^x^i^O;^,^^) 

p=0 
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by the g-binomial theorem valid for \x\ < 1 and rearranging the absolutely convergent sum 
using I = p — k. From this we obtain the recurrence relation (3.12). The general estimate 

\2<Pi{q p ,0;g;g,xg p )| < 



and Gronwall's inequality (3.9) prove that sum is absolutely convergent around 0. □ 

Choose r > so that j u -i(x; q) has no non-zero zeros in the ball with radius r and the 
origin as centre, which is possible since i<^i(0; qx 2 ; q, q u x 2 ) equals 1 at x = and defines 
an analytic function for \x\ < q~ x l 2 . Using lemma 3.3 we obtain for m G Z and n G Z + 

(3.13) — / z™ J ^4dz={°> ---°ddorm>-n-2, 
2m J\ z \ =r j v _i{z;q) { 1, m = -n-2. 

The coefficients dk of lemma 3.3 for n = can be interpreted as the moments of the 
moment functional £_ defined by 

1 f jJz~ l :q) f 0, m G Z + odd, 

(3.14) £_(x m ) = — <b z m Ju K . _\ q \ dz=\ + 

2m J\ z \=i ju-i{z 1 ;q) { d m/2 , m G Z + even. 

In §5 we consider the orthogonal polynomials for £_ from which some properties for j„(x; q) 
can be derived. 

Define the following strong moment functional C for v > on the space of Laurent 
polynomials by 

for any Laurent polynomial p(z) = Xl^Ln c v zV ■> n < m, n,m E T,. Note that C is indepen- 
dent of the choice of R, respectively r, as long as J v -\{x\ q), respectively j u -i(x; q), has no 
non-zero zeros in the ball with radius R, respectively r. All moments of £, both positive 
and negative, are well-defined due to lemmas 3.2 and 3.3. 

The moments of strong moment functional C and the moments of the moment func- 
tional C± defined in (3.11) and (3.14) are related by C + (x n ) = C(x n ), n G Z+, and by 
C-(x n ) = —C(x~ 2 ~ n ), n G Z_|_. 

Theorem 3.4. Let v > 0. The Laurent q-Lommel polynomials h UjU (x; q) defined by (1.7) 
are orthogonal Laurent polynomials with respect to the strong moment functional C, cf. 

(3.15) . Moreover, also the Laurent polynomials x~ 1 h ntl/ (x; q) are orthogonal with respect 
to C. Explicitly, 

C(h n ,„(x; q)h mt „(x; q)) = - » L{x~ x h n ^{x\ q)x~ 1 h m ^(x; q)) = -<5 n>m . 

Remark, (i) This result corresponds nicely with the fact that the Laurent g-Lommel poly- 
nomials correspond to a sequence of lacunary orthogonal Laurent polyomials, cf. (2.2), 
(2.3). 
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(ii) Since C(x 2 ) = -1 we see that £ is not a positive definite strong moment functional. 

Proof. The asymptotically well-behaved solutions J l/+n (x~ 1 ; q) and j u + n (x; q) of the recur- 
rence relation (1.7) are expressible in terms of the Laurent polynomials h njU (x; q) and the 
associated Laurent polynomials h n -i jU +i(x; q), cf. proposition 3.1. From this we obtain 
for any m E Z the expressions 

/o i r\ mJv+n{,% ] q) m Jv (^ 1 (?) , , \ m , / \ 

(3.16) x — — -=x — rh n , u (x;q) - x h n - lt „ +1 {x;q), 

J v -i{x L ;q) J v -i{x L ;q) 

(3.17) x - = x ?h njV {x;q)-x h n -i tV+1 {x;q). 

Since we obviously have 

— / z m h n - 1)U+1 (z]q)dz = — I z m h n - 1 ^ +1 (z;q)dz, 
2ttz J\ z \=^ 2tti J lzl=r 

we get from the combination of (3.16), (3.17), (3.10) and (3.13) the relations 

0, — n < m < n, 

{q v ;q)nln m = n, 

0, — n < m < n, 
— 1, m = —n — 1, 



C(x m h n Ax;q)) = 
C^x^K^x-q)) 



This proves the orthogonality. 

It remains to calculate the norm. From (3.1), (3.2) we see that the coefficient of x n in 
h n ,v{x',q) equals (q u ; q) n and that the coefficient of x~ n ~ 1 in x~ 1 h rtjU (x;q) equals 1. □ 

4. Orthogonal q-Lommel polynomials associated with the positive moments 

In this section we consider the orthogonal polynomials for the moment functional C+, 
cf. (3.11), which corresponds to the positive moments of the strong moment functional C. 
These polynomials are (/-analogues of the Lommel polynomials h njll (z), cf. (1.4). 

We consider the following three-term recurrence relation, 

(4.1) p n+1 (x) = x(l - q v+n )p n {x) - XnPn-^x), A 2n = q n , A 2n+1 = q" +3n+1 

with initial conditions p-i(x) = 0, po(x) = 1. Note that we can write the recurrence 
coefficient A n in closed form as g( i/ + n )(L( n + 1 V 2 _l - L n / 2 J)+|_ n / 2 J ^ w \ iere |^ a j denotes the greatest 
integer less than or equal to a G HL So the recurrence relation (4.1) depends on whether 
n is odd or even. Favard's theorem implies that these polynomials are orthogonal with 
respect to a positive definite moment functional for v > 0. Taking q f 1 in (4.1) after 
replacing x by 2z/(l — q) yields the three-term recurrence relation (1.4) for the Lommel 
polynomials, so that we have (/-analogues of the Lommel polynomials. The recurrence 
relation (4.1) has been guessed using the explicit form for the positive moments of £, i.e. 
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the moments of C + , obtainable from lemma 3.2 and calculating the first few terms of the 
recurrence relation (4.1) using Mathematica. 

The monic orthogonal polynomials satisfy a recurrence relation of the type 

r n+1 (x) = xr n (x) - fi n r n -i(x), r_i(x) = 0, r (x) = 1, 

with fi n > for all n G N and Yl^Li A*n < °°- This type of orthogonal polynomials has 
been studied by Dickinson, Pollak and Wannier [10], by Goldberg [13], who corrected some 
of the results of [10], and, from the point of view of continued fractions, by Schwartz [19]. 
See also Chihara [7, Ch. IV, thm. 3.5]. The support of the corresponding orthogonality 
measure, which is uniquely determined, is a purely discrete denumerable bounded set with 
only one accumulation point at zero. This result can also be obtained by remarking that 
the Jacobi matrix J for the corresponding orthonormal polynomials defines a self-adjoint 
operator J:£ 2 (Z + ) — > £ 2 (Z + ), which is an operator of trace class. Since the spectral 
measure of J is the orthogonality measure for the orthogonal polynomials r n , the result 
follows from standard facts on the spectral measure of a self-adjoint trace class operator. 
Moreover, for the orthogonal polynomials in this class we have the asymptotic behaviour 
of the form lin^^oo x~ n r n (x) = f(x) for an analytic function / in C\{0}, cf. [10], [13], 
[19]. 

By pn we denote the associated orthogonal polynomials, i.e. the polynomials satisfying 

(4.2) p^\x) = x(l-q^ +n )p^_ 1 (x)-X n p^_ 2 (x), \2n = q n , A 2n+1 = q» +Sn+1 

with initial conditions p^\(x) = 0, Pq (x) = 1. 

The following proposition is a (/-analogue of the identity (1.3) relating the Bessel func- 
tions and Lommel polynomials. 

Proposition 4.1. For n G Z + the polynomials defined by (4.1) and (4.2) satisfy 

Pn{\) Ju Or; q) ~ p£ ! {\) J v -i Or; q) = ? * L<*+D/*J (»+") J v+n ( xq h L(«+D/2j . q) 

where J u (x;q) denotes the Hahn-Exton q-Bessel function (3.3). 

Proof. The left hand side is a solution of the three-term recurrence relation 

1 - q v+n 

(4.3) a n+ i = a n — A n a n _i. 

x 

The right hand side satisfies the same recurrence relation (4.3). To see this we use for even 
n the relation 

(4.4) l^Lj^q) - Ju-^q) = q^ v+1 ^J v+1 (xq^;q), 

x 

and for odd n we use the relation 

(4.5) Ju(x;q) - q^' 1 ^ J v _ x {xq~^ ; q) = J u+1 (x;q). 

x 
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These identities can be checked straightforwardly by comparing the coefficients of the 
powers of x on both sides of (4.4) and (4.5). 

Since p n (^ _1 ) and p^l 1 (x~ 1 ) are linearly independent solutions of (4.3) we obtain the 
proposition after checking the equality for n = 0, which is trivial, and for n = 1, which is 
(4.4). □ 

The polynomials defined by (4.1) turn out to be the orthogonal polynomials with respect 
to the moment functional £+. For more information concerning the zeros of the Hahn- 
Exton g-Bessel function, which play a role in the following theorem, the reader is referred 
to [17, §3]. 

Theorem 4.2. We have the following orthogonality relations for v > for the polynomials 
defined by (4.1); 



k ' ^- lJ «r 1 ) 2 ->uor 1 ;9) 

? (n+^)L(n+l)/2j 
+Pn(0) Pm (0) = 5 n , m 1 _ y . 

Here are the positive simple zeros of the Hahn-Exton q-Bessel function J iy -i(x;q) 

numbered increasingly. All weights are positive. 

Proof. We start by establishing a complex orthogonality similarly to the previous section 
following Dickinson's method [9]. For this we need the expansion 

(4 6) q h\(n+l)/2\(n+v) Jv+n{Xq^ \q) = \ (n+1)/2 \(n+v) X V CfcX 2fc 

J„-i(x;q) (q»;q) n+1 ^ Q 

which is absolutely convergent for small x. Moreover, cq = 1. This is proved similarly as 
in lemmas 3.2 and 3.3. 

Let R > be smaller than the smallest positive zero j\~ x of J v -\{x\ q), then we obtain 
from proposition 4.1 and (4.6) for < m < n 

I Z m Pn {z) T MZ ~ 1 \ q \ dz= I ^ g H(n+D/2j (n+u) Jv+niz^qV L( " +1)/2J j g) ^ 

J\z\=ji Ju-i{z 1 ;q) J ]z]= i_ J„-i{z ^q) 

0, < m < n, 

2mq^ +n ^ n+i y 2 i(q v ;q)-l 1 , m = n, 

since §, , ± z m p^_ 1 {z)dz = 0. The leading coefficient of p n is {q v ; q) n , as can be read off 
from (4.1), and so we get the complex orthogonality relations 

(4.7) £+(p n p m ) = — <b Pm(z)p n (z)- p— ^ — rdz = S r , 



2 m J M= ^>^>J v _ 1 {z-\qr ^ (l-q»+ n ) ' 
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The considerations given at the beginning of this section show that we can rewrite (4.7) 
as a sum over the zeros of the Hahn-Exton g-Bessel function J v _\{z\ q) and possibly zero. 
The residues at the pole (j^ -1 ) -1 of the left hand side of (4.7) equals 

To see this we note that ^-iC/fc -1 ; g) 7^ since the zeros of J v -i(x;q) are simple, cf. 
[17, lemma 3.3], and that J v {j v k ~ x \ q) ^ by the interlacing property of the zeros of the 
Hahn-Exton g-Bessel function [17, thm. 3.7]. The positivity of the corresponding mass 
follows from the fact that Jv{j k ~ > 0) an d J' v -iUk~ ! •?) have opposite signs, which follows 
from the Fourier-Bessel orthogonality relations for the Hahn-Exton g-Bessel function [17, 
prop. 3.6], or from the fact that the zeros of the Hahn-Exton g-Bessel functions J l/ (x;q) 
and J u+ \{x; q) are interlaced as described in [17, thm. 3.7]. The mass at — (j^ -1 ) -1 yields 
the same weight. 

The set of mass points k e N, has zero as the only point of accumulation, 

so that zero may occur as a mass point as well. This happens if J2T=o \Pk(^)\ 2 < °°? 
where p n are the corresponding orthonormal polynomials, cf. e.g. [5, thm. 2.8]. Now the 
orthonormal polynomials p n are given by 

(4-8) p n (x) = [^Y^r) q-^ n+1)/2i{n+ ^Pn(x). 

Moreover, M(p n p m ) = 

^n,mi where A4 is the moment functional given by 

M(p) = 1 -^f I p(z) ^ {Z '\ q \ dz = (1 - q»)£ + (p). 

From (4.1) with j = 0we see that p2 n +i(0) = and that £>2n(0) satisfies a simple two- 
term recurrence relation from which we get P2n(0) = (— l) n g n ^ +1 )+2 n ( n_1 ). Combining 
this with (4.8) shows that p 2 n+i(0) = and 

i _ n u+2n \ 1/2 



fen(O) = (-1)" ( 1 _ qU ) q~'- 



^ OO 

* * ■ "+2n\ nv-\-n(n— 1) 



Hence, 

oo 

P = £fc(o)l 2 = ^- 7 £(i-<f +2 ") 9 

fc=0 q n=0 

and this sum is an absolutely convergent telescoping series so that p = (1 — q v )~ 1 - Con- 
sequently, M. has a mass point at zero with weight p~ x and C+ has a mass point at zero 
with weight 1. □ 

From the explicit orthogonality relations of theorem 4.2 we see that the orthogonality 
measure for p n {x) is supported in [— l/ii _1 , \j j\~ x \. On the other hand, from the explicit 
values of the recurrence coefficients for the orthonormal polynomials p n , which are easily 
obtained from (4.1) and (4.8), and the bound on the spectrum from [21, (1.3) with n — > oo], 
which is Gershgorin's theorem for the Jacobi matrix, we see that the orthogonality measure 
is supported in [-N, N] with N < 2/(1 — q u ). So we obtain the following corollary after 
shifting v by 1. 
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Corollary 4.3. Forv > — 1 the first positive zero j± ofJ l ,(x;q) satisfies j\ > (1 — q l/+1 )/2. 

5. Orthogonal polynomials associated with the negative moments 

In this section we consider the orthogonal polynomials for the moment functional £_ 
related to the negative moments of the strong moment functional C introduced in (3.14). 
In subsection 5.1 we introduce the three-term recurrence relation for the polynomials which 
we study. The three-term recurrence relation has been obtained by calculating the first few 
recurrence coefficients using lemma 3.3 with n = using Mathematica and then guessing 
the general result. In subsection 5.1 we give explicit expressions for these orthogonal 
polynomials and the associated orthogonal polynomials in terms of Al-Salam-Chihara 
polynomials. From the explicit expressions we can determine the asymptotic behaviour 
of the (associated) polynomials as the degree tends to infinity in terms of the function 
j v {x; q). In particular, we obtain the Stieltjes transform of the orthogonality measure. In 
subsection 5.2 we use the Stieltjes transform to obtain information on the zeros of j v {x; q) 
in a similar way as in Ismail [15, §4], see also [2, §4], and to give explicit orthogonality 
relations. In subsection 5.3 we give a different derivation of some of these results in the 
special case v = 1/2, which turns out to be related to known orthogonal polynomials [2], 
[20] . Comparison of these two approaches yields a summation formula for a one-parameter 
terminating 3</? 2 -series. A special case of this summation formula is the evaluation of the 
continuous g-Hermite polynomials at a special point outside the spectrum. 

5.1. Explicit expressions for orthogonal polynomials. We investigate the monic 
orthogonal polynomials satisfying the three-term recurrence relation 

(5.1) P n+1 {x) = xP n {x)-\ n P n - 1 {x), \ 2n = q n , X 2 n+i = q n+ ^ 

with initial conditions P-i(x) = 0, Po(x) = 1. By Favard's theorem the polynomials P n 
are orthogonal with respect to a positive definite moment functional for i/Gl Moreover, 
the polynomials P n fit into the same class of Dickinson, Pollak and Wannier [10], Goldberg 
[13] and Schwartz [19] described at the beginning of the previous section. 

The polynomials P n are even functions of x for even n, and odd functions of x for odd 
n. Introduce 

P 2n (x) = R n (x 2 ), P 2n+1 (x) = xS n (x 2 ), 

so that the monic polynomials R n and S n satisfy the three-term recurrence relations (see, 
e.g., Chihara [7, p. 45]) 

Rn+l(x) = (x — \2n — Mn+l)Rn{x) ~ Mn-l MnRn-1 (x) , 
S n +l(x) = [X — A 2n +1 — ^2n+2)S n {x) — \2n^2n+lS n -l(x), 

with initial conditions Ro(x) = 1, Ri(x) = x — q v and S-i(x) = 0, So(x) = 1. A simple 
computation from (5.1) gives the recurrence coefficients for the polynomials R n 

A2n + A 2n +1 — \ ,, . e n ) ^2nA2n-l ~ Q , U >\i. 

y q if n = 0. 
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For the recurrence coefficients of S n we find similarly 

A 2 n+i + A 2 n+2 = (q + q u )q n , A 2 nA 2 n+i = q 2n+ " n>0. 

The recurrence coefficients of R n and S n decrease exponentially. 

Consider the monic polynomials u n (x; a, b; q) satisfying the recurrence relation 

(5.2) u n+ i(x; a, b; q) = (x - aq n )u n (x; a, b; q) - b 2 q 2n ~ 2 u n -i(x; a, b; q), 

U-i(x) = 0, uo(x) = 1, which are studied in [20], then S n (x) = u n (x;q + q u ,q^ +2 ^ 2 ;q). 
For R n we have to be a little bit more careful, since for n = one of the recurrence 
coefficients behaves differently. However, R n is still a solution of the recurrence relation 

(5.2) with a = 1 + q v and b 2 = q u+1 , but it satisfies the different initial condition R\{x) = 
x — q v = u\{x) + 1. Such polynomials are known as co-recursive polynomials [6] and can 
be expressed as 

R n (x) = u n (x; 1 + q\ g^+D/ 2 ; q) + ^(j; 1 + q\ q^' 2 - q). 

The associated polynomials corresponding to the recurrence relation (5.2) are given by 
(x; a, b; q) = u n (x; aq, bq; q) = q n u n (x/q; a, b; q), so that 

R n (x) = u n {x- 1 + q\ q^' 2 - q) + q n ~' u n _ x {x , 'q; 1 + q\ q^/ 2 ; q). 

An explicit expression of the polynomials u n (x; a, b; q) in terms of Al-Salam-Chihara 
polynomials is given by Van Assche [20, thm. 2]; 

u n (x; a, b;q) = J2 — r~\ Pk(~a; q; -aq n ~ k+ \ b 2 q 2 ^ + \b 2 /q). 

fc=o [q ' q)k 

Here P n (x;q;a,b 7 c) are Al-Salam-Chihara polynomials, cf. [1, (6.1)], which satisfy the 
recurrence relation 

(5.3) P n+1 (x;q;a, 6, c) = (x - aq n )P n (x; q; a, 6, c) - (c - 6<? n_1 )(l - q n )P n -i(x; q; a,b, c). 

More information, including the orthogonality relations, concerning the Al-Salam-Chihara 
polynomials can be found in Askey and Ismail [5, §3]. 
So we obtain the explicit expressions 

71 T n-k fc(fc-l)/2 

(5.4) s n (x) = ^-r—, p k(-(Q + Q u );q; -(i + ^"V - ** 2 . q 2 ^- k ^+ 3 , <f +1 ) 

fc=0 

and 

n T n-k fc(fc-l)/2 

^0*0 = £ — ^— P fc (-(i + O; «;-(i + o«f _ * + V (B_ * )+,/+2 , g ") 

fc=0 W,9Jfc 

n-l-fc fc(fc+l)/2 

+ £ jt^t P*(-(l + «"); <z; -(1 + ^)? n " fc , ^- 1 - fc )+^+ 2 , <f ) 



fc=0 



r n — k n k(k — l)/2 



= x n + Yl X (l q)k " ^(-(! + OS «; "(I + ^)9 n - fc+1 , q 2 ^ + » +2 , q v ) 



k=i 



+ (i - ? fc ) p fc _i(-(i + O; «; -(i + 9")« n -* +1 , 9 2 < n - fc > + " +2 , o] 



ORTHOGONAL POLYNOMIALS AND LAURENT POLYNOMIALS 



17 



A generating function for the Al-Salam-Chihara polynomials is, cf. [1, p. 23], 



Mx 2] = f P (x . Q . a b c) ^_ = (^;g)oo(^;g) 



oc 
oo 



n=0 

where 1-az + bz 2 = (l — az)(l — {3z) and 1-xz + cz 2 = (l—'yz)(l — 6z). Take x = —(l + q v ) 
and c = q v so that 7 = — 1 and 5 = — q v . Consequently (1 + z)$>(z, — (1 + q u )) is the 
generating function for x = —(q + q v ) and c = q u+1 . Hence, 

P n (-(1 + q»);q; a, 6, g") + (1 - q n )P n -i (-(1 + g"); q; a, 6, g") = P n (-(q + q v ); q; a, 6, 
and thus 

" n-fc„fe(fe-l)/2 

(5.5) i? n (x) = £ 9 : + 9"); 9; -(1 + «")« n -* +1 , ff 2 <»- fc > + "+ 2 , 

fc=0 W,9jfc 

Now that we have the explicit expression for the polynomials P n defined in (5.1) at 
hand, we can determine the asymptotic behaviour, which is related to the function j u (x; q) 
introduced in (3.5). 

Proposition 5.1. For the orthogonal polynomials P n (x) defined by (5.1) we have for every 
x E C 

lim x n P n (l/x) = x 1 - / j v . 1 (x;q). 

n — >oo 

Proof. We follow the proof of theorem 2 of [20] . For this we need the continuous g-Hermite 
polynomials H n {x \ q) introduced by Rogers in 1894. The three-term recurrence relation 
is 

(5.6) H n+1 (x I q) = 2xH n (x \ q) - (1 - g n )# n _i(z I ?) 

with initial conditions H-\(x \ q) = 0, H (x \ q) = 1, cf. Askey and Ismail [4, §6]. From 
(5.3) and (5.6) we obtain, cf. [20, thm. 2], 

(5.7) Um ^-vr^w^.m - (-i)^-W^ I q). 

n^oo \ zo / 

Using this limit relation and dominated convergence we obtain 

lim x n R n {l/x) = lim x n S n (l/x) 

n — >oo n — >oo 

and hence 

00 ( ->\k k(k — l)/2 1 

(5.8) lim x n P n (l/x) = y [ ~ L) , q , x 2k qi^H k (Uqh^) +q h(^)) \ q \ 

n^oo ^ (g;g)fc ' V2 / 
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To see that the right hand side of (5.8) equals x 1 u j v -i{x;q) we insert the explicit 
expression, cf. [4, (6.1), (3.1)], 



H k (±(x + x- 1 )\q)=J2 

1=0 



(T,q)k k -2i 



(q;q)i(q;q)k- 



for (5.8). Interchanging summations and introducing m = k — l shows that 

(5.8) equals 



( i\L|l(l+l) T 21 °° ( 1\m im(m-l) 
\ L ) q 2 X q 2 x 2m q m(l+is) 



1=0 vi' i/' m=0 

oo 



r-i v«I / ( /+1 )t- 2/ 

by use of [11, (1.3.16)]. □ 

Observe that the continuous g-Hermite polynomials are orthogonal on the interval 
[— 1, 1], so that the inequality 2 < q2^- v ) + g 2 !) shows that the variable of the continu- 
ous g-Hermite polynomial in (5.8) lies outside the support of the orthogonality measure for 
the continuous g-Hermite polynomials, except when v = 1 in which case it is an endpoint 
of the interval. 

The Stieltjes transform of the orthogonality measure \x for the orthogonal polynomials 
P n can be obtained from 



(5.9) / = lim 



Jr z ~ 



t n^oo P n (z) 



where Pn are the associated polynomials, cf. Askey and Ismail [5, thm. 2.4] and further 
references therein. 

So let us now consider the associated monic polynomials P n ^ satisfying 

(5.10) P^Ux)=xP£\x)- ln p£\(x), P7(*) = 0, P (1) (*) = l, 

where 7™ = A n +i is defined in (5.1). These polynomials can be determined in a similar 
fashion as before. Because of the parity of these polynomials, we again set 

and the monic polynomials T n and U n then satisfy the recurrence relations 
T n+1 (x) = (x - 7 2n - <y 2n +i)T n (x) - n/2n-i1/2nT n -i(x), 

U n+1 (x) = (X - 7 2n+ l - ^2n+2)U n {x) - 7 2 n72n+l U n - 1 (x) , 
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Tq(x) = 1, Ti(x) = x — q and U-i(x) = 0, Uq{x) = 1, where 

/ (q + qnq n ifn>0, 2n+I , 

72n + 72n+l = < ., n , I2nl2n-1 = q , ™ > 0, 

^ g it n = U. 

and 

72n+l + 72n+2 = (1 + 72n72n+l = <? n+ " +1 , U > 0. 

Hence 

t/ n (x)= Wn (x; ? (l + ^),^+ 3 )/ 2 ;Q) 

n — A; k(k— 1)/2 

= E ^(-«(! + «"); 9; -(i + rt n - fc+2 , <z 2 <"- fc > + " +4 , <f +2 ). 

fc=0 

The polynomials T n are again co-recursive polynomials for the recurrence relation (5.2) 
with a = q + q v and b = g( I/ + 2 )/ 2 , with T±(x) = ui(x) + q v ', and thus 

T n (x) = g + q\ qW*; q) + q'+^u^x/q; q + q\ ^+ 2 V 2 ; q). 

From the generating function of the Al-Salam-Chihara polynomials we find 

P n (-q(l + q v );q;a,b,g u+2 ) = 

P n (-(q + «");«; «, &, + - + O; «; «, 6, <z" +1 ), 

so that 

n n-k k(k-l)/2 

The proof of the following proposition is analogous to the proof of proposition 5.1. 
Proposition 5.2. For every x E C we have 

lim x n P^ l \l/x) =x-"j v (x;q). 

n^oo 

5.2. Zeros of j v (x;q) and orthogonality relations. Combining propositions 5.1 and 
5.2 and (5.9) shows that the Stieltjes transform of the orthogonality measure ji for the 
polynomials is 



(5.11) 



dfi(t) ju(l/z;q) 



z-t j v _ 1 {l/z;q) 



for all z £ supp(d^), cf. [5, thm. 2.4]. From the Stieltjes transform we can derive the 
orthogonality relations for the orthogonal polynomials P n defined in (5.1). We start with 
an investigation of the zeros of j v {x; q). It turns out that the zeros of the function j u (x; q) 
behave similarly to the zeros of the (Hahn-Exton g-)Bessel function for v > — 1. The 
method of proof largely follows Ismail's investigation [15] of the roots of the Jackson q- 
Bessel function, see also [2, §4]. 
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Theorem 5.3. Let v e ffi and let the function j v {x;q) be defined by (3.5). 

(i) The functions j v {x; q) and j u+ i(x; q) have no common zeros, except possibly x = 0. 

(ii) The zeros of x~ u j v (x; q) are real, simple and symmetric with respect to x = 0. There 
are infinitely many of them and their only point of accumulation is oo. 

(iii) The zeros of x~ v f v {x;q) and x~ l/ ~ 1 j^+i(x; q) interlace. Moreover, the smallest positive 
zero of x~ u j v (x;q) is smaller than the smallest positive zero of x~ u ~ 1 j u+ i(x;q). 

Proof. First we prove (i) by use of an equality for i^i-series. The relation 

— % 

(5.12) i¥>i(0; c; q, z) - iv?i(0; c; q, qz) = j— - iv?i(0; cq; q, qz) 

can be proved directly or it can be obtained from one of Heine's contiguous relations for 
the 2V 7 i- ser i es ? cr *- [11, ex. 1.9(iv)]. Take c = qx 2 and z = q v+1 x 2 in (5.12) to get from 
(3.5) 

(5.13) j v (x; q) - x~ 1 j u+1 (x; q) = -q 1+u/2 x 2 j iy (x^/q; q). 
Substituting c = q u+2 x 2 , z = qx 2 in (5.12) and using (3.5) gives 

(5.14) 3v+i( x '>Q) -Q~*"xju(xy/q;q) = -q^ {l - u) x 2 j v+1 {x^q). 

If 7^ a is a zero of j u (x;q) and j u +i(x;q), then (5.13) implies that a^/q is a zero of 
j v {x; q). Next (5.14) implies that a^/q is a zero of j u+ i(x; q) as well. So aq k l 2 , k G Z + , are 
zeros of the analytic function x~ u j l/ (x;q), which implies that this function is zero. This 
contradiction proves (i). 

To prove (ii) we recal that the orthogonality measure d/i is supported on a bounded 
denumerable discrete set with zero as the only point of accumulation. So let d\i have mass 
Ak at the points {tk}kLn then (5.11) is 



oo 



(5-15) EA-= ^/u^V " 6t - 

The zeros of x 1 ~ u j iy -i(l/x; q) correspond precisely to the non-zero poles tk of the left 
hand side. So the zeros are real and simple. Since {tk}%%i has zero as the only point of 
accumulation, the only point of accumulation of the zeros of f v -\{x\ q) is infinity. 
To prove (iii) we consider the (positive) mass of d\i at a non-zero t^, 

0<A k = -t 2 Ml/tk > q) 



k 3l-i(Vt k ;qy 

So jv(a;q) and j' u _ 1 (a;q) have opposite signs for ^ a a zero of j v -i(x;q). If < 
a < b are two consecutive zeros of j v -i(x;q), then j' iy -i(a;q)j' u _ 1 (b;q) < 0. Hence also 
jv(a;q)j u (b;q) < and j v (x;q) has at least one zero in (a, b). In the interval (1/6, 1/a) 
both sides of (5.15) are differentiable, and the derivative of the left hand side is strictly 
negative. If j u (l/z; q) has more than one zero in (1/6, 1/a), then the derivative has a zero 
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in that interval. Thus j u (x; q) has precisely one zero in (a, b). This proves the interlacing 
property. 

Denote by x v k the positive zeros of jv(x; q) numbered increasingly; 

< x\ < x v 2 < . . . < x^ < x v j+1 < ... . 

Then it remains to prove that x\~ x < x\. Since x u j v (x; q) equals 1 for x = we get that 
j'u-ii^i -1 ; q) < and thus jv(x^~ ;q) > 0. So jv(x;q) has an even number of zeros in 
(0, x\~ l ), and the same argument as in the previous paragraph shows that this number is 
zero. □ 

The following proposition is the analogue of proposition 4.1 for the orthogonal polyno- 
mials P n and the functions j l> (x;q). 

Proposition 5.4. For n £ Z + the polynomials P n and P^ defined by (5.1) and (5.10) 

satisfy 

D /l w , , x / q m{m+hv) x 2m 3v{xq^-q), n = 2m, 

x x [ q m™+ 2 [.v ^'x^^-xixq^-.q), n = 2m-l, 

where j v {x; q) is defined in (3.5). 

Proof. It suffices to show that the right hand side satisfies (5.1) with x replaced by 
since the left hand side satisfies this equation and the cases n = (trivial) and n = 1 
(from (5.13)) are easily proved. For n = 2m this follows from (5.13) with x, v replaced by 
xq m l 2 , v — 1, and for n = 2m — 1 this follows from (5.14) with x, v replaced by xq^ m ~ 1 ^ 2 , 
v-\. □ 

In the proof of theorem 5.3 we obtained information on the orthogonality measure for 
the polynomials P n defined in (5.1). In the next theorem we describe the full orthogonality 
relations. This theorem can also be proved from proposition 5.4 analogously to the proof 
of theorem 4.2 from proposition 4.1. 

Theorem 5.5. Let v e M, denote by , k £ N, the positive zeros of the function 
j u -i(x;q) defined in (3.5). Then for the polynomials P n defined by (5.1) we have the 
orthogonality relations 



+ (i- q »- i )p n (o)p m (o) = 6 n , m] w+m+v) 



q¥( l +»\ n = 2l, 

q W+m+v)^ n = 2 Z + l, 



where the mass at x = only occurs for v > 1. All weights are positive. 

Proof. The only statements to be proved concern the norm and the weight at x 
Denote the squared norm of P n by || P n || 2 , then (5.1) implies, cf. [10, (7)], 

II Pn || 2= A n || P n -1 || 2 ==>■ || P n || 2 = A n . . . Ai || 1 || 2 . 
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Together with the explicit value for A n in (5.1) the statement on the norm follows if we 
prove || 1 || 2 = 1. The value of || 1 || 2 can be read off from the Stieltjes transform (5.11) as 
the coefficient of z~ x on the right hand side and lemma 3.3 for n = and x = z~ x shows 
that it equals 1. 

The weight at x = equals p, where p~ l = Y^=o -fn(O) 2 and P n denote the orthonor- 
mal polynomials, cf. [5, thm. 2.8]. From (5.1) we compute P 2 n+i(0) = 0, -P2n(0) = 
(— i) n qi n ^+"("- 1 )/ 2 so that for the orthonormal polynomials we have 

p 2n(0) = yo) = ( _ 1)v >(,-i)/ 2) 

V A1A2 • • • A2n 



and thus 



OO OO s 

E^ 2 ( ) = E^" 1) = 



n=0 n=0 



OO if V < 1, 

(1-q"- 1 )- 1 ifi/>l, 



so that there is a mass 1 — q v at the origin whenever v > 1. □ 

Again, as in the proof of corollary 4.3, using [21, (1.3) with n — > 00] shows that the 
orthogonality measure for the _P n is contained in [—AT, AT] with iV < 1 + q v l 2 . Shifting v 
to v + 1 we get the following corollary. 

Corollary 5.6. The first positive zero x\ ofj 1/ (x;q) satisfies x\ > (l + g(^+ 1 )/ 2 ) . 

5.3. The case v = 1/2. In the simple case z/ = 1/2 we have \ n = q n l 2 . For simplicity 
we take p = q 1 ! 2 so that the recurrence relation (5.1) can be rewritten as 

(5.16) P n+1 (x)=xP n (x)-p n P n - 1 (x). 

We consider the generating function G(z, x) = J2^=o Pn(x)z n . Multiply (5.16) by z n+1 
and add all the terms from n = to infinity, then we get 

1 z 2 p 

G(z : x) — 1 = xzGiz, x) — z 2 pG(pz, x) G(z, x) = G(zp, x). 

1 — xz 1 — xz 

Solving the p-difference equation with respect to the condition G(0, x) = 1 gives by itera- 
tion 



We use the p-binomial theorem, 



(zx;p) k+1 ^ (p;p)n 
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in (5.17). Changing the summation index n to j — 2k gives 

3 fc+1 ; 

(p;p)j-2k 



( h-\-\ \ 

G(z,x) V V ' P)J -" 



k=0j=2k 
oo b'/2j 



fc j_ 2fc fc a (p fc+1 ;p)j- 2 fc 



V:' V(-1)V -p" 
Next identify the coefficient of 2 n and use (p fc+1 ; p)j-2k = (PiP)j-k/(p'iP)k to find 
(5-18) P„(x) = V (-l)*x"- 2 V a r ( ff )n , fc • 

£z£ (p;p)k{p;p)n-2k 

These polynomials are a special case of orthogonal polynomials associated with the 
Rogers-Ramanujan continued fraction; they correspond to the case a = 0, b = p and q = p 
in Al-Salam and Ismail [2] and (5.18) corresponds to [2, (3.7)]. These polynomials are 
also the special case u n (x) in Van Assche [20] with a = 0, b = q and q 2 = p, and (5.18) 
corresponds to [20, (2.7)] after observing that for the Al-Salam-Chihara polynomials in 

(5.3) we have 

P 2n+1 (0; q- 0, 6, c) = 0, P 2n (0; q; 0, 6, c) = (-l) n c n (-; Q 2 ) n ( ? ; q 2 ) n . 

c 

The associated polynomials P^ satisfy the recurrence relation 

(5.19) Pi%(x) = xP£\x) -p^P&ix), 

with P^ = and Pq 1 \x) = 1. Replace x by x/y/p in (5.19), then the monic polynomials 

p n / 2 P n (x/ y/p) satisfy the recurrence relation (5.19) so that Pn\x) = p n / 2 P n (x/ ^/p). 
In case v = 1/2 we have two different expressions for the same polynomials. From 

(5.4) and (5.18) we obtain the following summation formula for the Al-Salam-Chihara 
polynomials, < k < n, 

(5.20) JV-( 9+9 *);,;-(l +9 - *),-'«, ff*-^*,,*) = (^fe^gWlzj . 

(9 2 ;9 2 )fc(? 2 ;? 2 )2n+l-2fc 

The Al-Salam-Chihara polynomials are expressible in a 3^2-series as proved by Askey 
and Ismail [5, §3.8]. Explicitly, the following connection between the original notation of 
Al-Salam and Chihara [1] and the notation of Askey and Ismail [5] holds; 



a 



-k 



(q;q) 



-P k (2ax;q; (7 + 5)a^5a ,a ) = S k (x;^,S | q) 



k 



(5.21) 



7 3<^2 , _ ; q, q 



(q;q)k \ n/S, o 

where x = (y + y~ 1 )/2. 
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Corollary 5.7. The summation formula 

I 2 \ ( Q. k 7 C Q 1//2 ,C \ -\l1\k i 1/2 l/2\ / l/2\ 

(c ;q) k 3^2 I c2 ;q,qj ={cq 1 ) (-g 7 ;g 7 )fc(c;g 7 ) fc 
Zio/ds /or G Z + and c G C. 

3 11 

Proof. In (5.20) we use (5.21) with the parameters a = — g 1 , x = (g 3 4- q~i)/2, 7 = 
q n ~ k +4 ? 5 = g n_fc +i to get the result of the proposition for c = Replace n — k by 

m in this result, so that we have proved the corollary for c = g m+1 , m G Z+. Since both 
sides are polynomial in c, the result follows for arbitrary values of c. □ 

Remark. Comparison of (5.5) with (5.18) instead of (5.4) with (5.18) leads to the same 
corollary. The same result is also obtained if we work out the different expressions for the 
associated polynomials in case v = 1/2. 

Proposition 5.8. Consider the monic orthogonal polynomials given by (5.16) and the 
associated polynomials given by (5.19). Then for every x G C we have 

lim x n P n (l/x) = F(x), lim x n P£ ) {l/x) = F(xJp) 

n^oo n^oo 

where 

F K X ) = 2^ — 7 \ =0^1 ;p, -x p . 

Proof. Straightforward by letting n — > 00 in (5.18) after changing x to 1/x and multipli- 
cation by x n . □ 

From proposition 5.8 and propositions 5.1 and 5.2 for v = 1/2 we obtain the equalities 

(5.22) o¥>i(-;0;g5, -x 2 g*) = x^j_ 1/2 (x; q), o<pi(-; 0; g^ , -x 2 q) = x~^j 1/2 (x;q), 

which gives two transformations of a o¥ : 'i" ser i es °f base g 1 / 2 = p in terms of i</?i-series 
of base g. In two special cases the left hand sides of (5.22) can be summed by the 
Rogers-Ramanujan identities, cf. e.g. [11, (2.7.3), (2.7.4)] and this gives explicit values for 
j-i/2( x 'i q) f° r x = =t*) Hq 1 ^ 4 and for ji/2(x; q) for a; = ±z, ±ig -1 / 4 . 

On comparing powers of a; in the first equality of (5.22) and using the expansion (5.8) 
with v = 1/2 for x j v -\{x; q) we find for the continuous g-Hermite polynomial the 
evaluation 

^Q(g^+g-)|g)=g- fc/4 (-g 1/2 ;g 1/2 ),. 

The same result is obtained if we compare powers of x in the second equality of (5.22). 
Moreover, this result is a limit case of corollary 5.7, cf. (5.7). 
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6. Orthogonality for the Laurent q-Lommel polynomials 

In this section we give a different form for the strong moment functional introduced 
in §3. The limit transitions (3.4) and (3.6) suggest to rewrite the strong moment func- 
tional C defined in (3.15) as a contour integral over the unit circle. This can be done if 
j v -\{l; q) ^ 0, since we have sufficient knowledge on the location of the zeros of j u -i(x; q), 
cf. theorem 5.3, and of J u -i(x;q), cf. [17, §3] and §4. A Wronskian type formula can be 
used to simplify the integrand. 

Lemma 6.1. Let r m {x), s m (x) be solutions of the recurrence relation (1.7), then the 
Wronskian r m (x)s m+ i(x) — s rn (x)r m+ i(x) is independent of m G Z. 

Proof. Multiply the recurrence formula for r m (x) by s m (x) and multiply the recurrence 
relation for s m (x) by r m (x). Subtract the resulting identities to find the result. □ 

Lemma 6.2. 

-i (^" 2 ;?)oo(^ 2 ;?)oo 



Ju(l/x;q)j v -i(x;q) - J v -i(l/x;q)j v (x;q) = x 



Proof J u+m (l/x; q)j l/+rn _ 1 (x; q) - J^ +m _i(l/x; q)j l/+m (x; q) is independent of m by pro- 
position 3.1 and lemma 6.1. Take m = to obtain the left hand side of the lemma and 
use (3.7) and m — > oo to see that it also equals 



, -l M x ,q)oo(qx 2 ;q) oo 

which proves the lemma. □ 

Lemma 6.2 implies that x(q; q)^ [j v (l/x; q)j u -i(x; q) — J u -\(l/x\ q)jv(x; q)] is a theta 
product, cf. Askey [3, §1]. 

Now we can rewrite the strong moment functional C with respect to which the Laurent 
g-Lommel polynomials are orthogonal, cf. theorem 3.4. 

Theorem 6.3. Let s > such that s is not a zero of J u -i(l/x;q) and j v -\{x\ q). For 
v > the strong moment functional £ defined in (3.15) equals 

C(n) = 1 I p f z \ (^~ 2 ;g)oo(^ 2 ;g)oo dz_ + 

2ni(q; q) oo J\ z \= s J„- 1 (l/z;q)j„- 1 (z;q) z 

where p is an arbitrary Laurent polynomial. Here , respectively x\~ x , denote the 
positive zeros of J u -i(x;q), respectively j v -i(x;q), numbered increasingly. N is defined by 
j 1 ^ 1 < s < jtf+i and N = 0, and so the sum over the zeros of J J/ -i(x;q) is empty, if 
f\~ x > s. M is defined by x 1 ^ 1 < s < x v ^^_ x and M = 0, and so the sum over the zeros of 
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j u -i(x;q) is empty, if x\ > s. The discrete weights in the first sum over k are positive 
and the discrete weights in the second sum over I are negative. 

Remark, (i) By choosing s = r, respectively s = 1/R, with r, R as in §3, we get M = 0, 
respectively N = 0. In §7 we show that for v sufficiently large we have N = M = for a 
suitable choice of s. 

(ii) The non-zero poles of the integrand in theorem 6.3 are simple. Indeed, if ^ a satisfies 
J u -i(l/a;q) = = j u -i(a;q), then lemma 6.2 implies that the numerator is zero as well. 
Moreover, a = q p / 2 for some p G Z, which is a simple zero of the numerator. There exist 
only finitely many of such values in the (possibly empty) interval [x"~ , l/j^ -1 ]. 

Proof. In the first contour integral in (3.15) we shift the contour integration from \z\ = 1/R 
to \z\ = s and in general we assume s < 1/R. We pick up residues at the simple poles 
z = , k = 1, . . . , AT, cf. theorem 4.2. For 1/R < s we are in the case N = 0. The 

second contour integral in (3.15) is shifted from |z|=rto|2;|=s. In general we assume 
r < s, otherwise we are in case M = 0. Here we pick up residues at the simple poles 
z = ±a^ _1 , 1 = 1,..., M. The residues are easily calculated. Next we take together the 
integrands of the contour integrals over \z\ = s using lemma 6.2 to prove the expression 
for C(p) in this case. The last statement follows from theorem 4.2 and theorem 5.5. □ 

Remark. The most natural choice for s in theorem 6.3 seems s = 1. This is motivated 
by the fact that on the unit circle there is a transition in the asymptotic behaviour of 
the Laurent g-Lommel polynomials, cf. (3.4), (3.6). Moreover, numerical experiments 
indicate that for m — > 00 the non-real zeros of the Laurent g-Lommel polynomials, cf. 
remark 2.1, are possibly dense on the unit circle. Of course, from (3.4), respectively (3.6), 
we see that the real zeros outside, respectively inside, the unit circle tend to the zeros of 
J 1/ -i(x~ 1 ; q), respectively j u -i(x; q). This corresponds precisely with the discrete set in 
the orthogonality measure of theorem 6.3 for s = 1. 

7. Laurent q-Lommel polynomials as 
perturbations of chebyshev polynomials 

Let us now return to the recurrence relation (1.7) which we rewrite as 

(7.1) h n +i, v (x; q) - (x' 1 + x)h n ^(x\ q) + h n -i, v (x; q) = -xq u+n h n ^{x; q). 

In this way, as q — > or as v — > 00 the Laurent polynomials h n ^(x; q) should be close to 
a solution of the three term recurrence relation 

(7.2) h n+1 (x; 0) - (x' 1 + x)h n (x; 0) + /i n _i(x; 0) = 0. 

The solution of this recurrence, with initial values ho(x;0) = 1 and h-i(x;0) = is given 
by h n (x;0) = (x n+1 — x~ n ~ 1 )/(x — x -1 ), which in terms of Chevbyshev polynomials of 
the second kind can be written as 

( x ~\~ x \ 
h n (x;0) = U n [ J, neZ + . 
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In this way the Laurent polynomials h nj „(x; q) can be considered as perturbations of the 
Chebyshev polynomials. We now do a perturbation analysis, much as is done for pertur- 
bations of orthogonal polynomials in [21]. In the spirit of the Liouville-Green approxi- 
mation (WKB method), we will consider (7.1) as a second order recurrence relation with 
non-homogeneous term —xq u+n h njU (x;q), even though this term depends on the desired 
solution h n ^(x; q). 

We solve this non-homogeneous recurrence relation by Green's method. We need 
the Green function Gi(n,m), which is the solution of the recurrence relation with non- 
homogeneous term 5 njm , i.e., 

(7.3) Gi(n + 1, m) — (x _1 + x)G\(n, m) + Gi(n — 1, m) = 5 n>m 
with boundary conditions 

(7.4) Gi(n, m) = 0, n>m. 

Clearly Gi(m, m) = Gi(m + 1, m) = and thus from (7.3) we find Gi(m — 1, m) = 1. For 
k > we find that ru{x) = Gi(m — k — l,m) is a solution of the homogeneous recurrence 
relation (7.2) with the same initial conditions ro(x) = 1 and r-±(x) = 0, hence 

G 1 {n,m) = U m - n - 1 (^-^- — V n < m. 



Now multiply (7.1) by Gi(n, m) and (7.3) by h n ^ v {x\ q) and subtract the obtained equations 
to find 

h n +i,v(x; q)Gi(n, m) - h njU (x; q)G 1 {n - 1, m) + h n ,„(x; q)5 n , m 

= h n ,„(x; q)G 1 {n + 1, m) - h n -i, u (x; q)G 1 (n, m) - xq" +n h njU (x; q)Gi(n, m). 

Add all the equations from n = to n = m and use the boundary conditions (7.4) to find 

m — l 

h , u {x; q)G 1 (-l, m) = h m ^(x; q) + x^2 q u+n G 1 (n, m)h n ^(x; q). 

n=0 

This gives 

(x -\- x~ 1 \ 771—1 ( x -V x~ 1 \ 
2 ) ~~ X ^ q u+n U m -n-i f J h n , v (x; q). 

From this relation we can deduce some useful properties. 
Lemma 7.1. Suppose x = e 10 with 9 G [0,2n), then 



(7.6) \h njU (x; q)\ < (n + 1) exp 



q v 



a-qy 
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and 
(7.7) 

For \x\ ^ 1 we have 



q v ( q v 

sin# h n ,v(x; q)\ < 1 + — exp 



2 / 2 q v \ 

\x n h n)U (x;q)\ < T exp -- , \x\ < 1, 

|1 — x z \ \\1 — x z \ 1 — q J 

\x~ n h n ^(x;q)\ < - — exp ( - — ^-^r-^- 
1 — x 2 \ 1 — x 2 1 — 



\x\ > 1. 



Proof. We use Gronwall's inequality, cf. (3.9); for non-negative A, c n , d n , (n > 0) we have 

n-l /n-l \ 

c n < d kCk =>- c n < A exp I ^ dk j • 

fc=0 \fc=0 / 

From the bound |L/ n (cos#)| < n + 1 and (7.5) we find 

n-l 

| h n , u (x; q)\<n + l + J2 Q" +k (n - k)\h kjU (x; q)\. 

k=0 

Hence taking c n = \h nj v {x\ q)\/(n + 1) in Gronwall's inequality gives 



'n-l 



\K,is(x;q)\ , ^ 

'—— < exp > (k + l)q 

n+1 \ r— ^ 



+k 



xk=0 



The desired inequality (7.6) then follows from YlT=o(^ + -0? = (1 — 2 - If* we use this 
inequality (7.6) and | sin 9 U n (cos9)\ < 1 in (7.5), then 



sin# h n ,„(x; q)\ < 1 + exp (^ (1 _ g)2 j X> + 



n-l 



which gives (7.7). The bounds away from the unit circle follow by using 



and 



\x n U n 



\x~ n U n 



X + x~ 



x n+l _ x —n—l 



X + X 



-1 



X 



X 



-1 



x — X 



x n+l _ x —n—l 



< 



< 



x — x~ 



1-X 2 V 



\l-x~ 2 V 



\x\ < 1, 



\x\ > 1, 



and by using Gronwall's inequality. □ 
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From these bounds we see that the Laurent Lommel polynomials have an exponentially 
increasing upper bound both inside the unit circle and outside the unit circle, and that on 
the unit circle the Laurent polynomials are bounded, except when x = ±1, in which case 
\h njV (x; q)\ = 0(n). This strongly suggests that in theorem 6.3 the choice s = 1 for the 
strong moment functional C is the most natural. 

The Laurent polynomial solution of (7.1) is not the only interesting solution. In §3 
we already obtained the minimal solutions j u + n (x; q) and J iy+n (x _1 ; q) on respectively the 
open unit disk and the exterior of the closed unit disk. The minimal solutions h~(x; 0) and 
h^(x; 0) of the recurrence relation (7.2) on respectively the open unit disk {z G C : \z\ < 1} 
and the exterior of the closed unit disk {z G C : \z\ > 1} are given by h~(x; 0) = x n and 
h+(x; 0) = x~ n . Our intention now is to find similar solutions u (x; q) satisfying 



lim h ( x '^) x ~ 

n— >oo ' 



on {z G C : \z\ < 1} and {z <E C : \z\ > 1} respectively. Such functions clearly exist, since 
by proposition 3.1 and (3.7) we see that 

(7.8) h+ v (x; q) = x v J l/+n (x- 1 ; q), h~ v (x; q) = ^J^^ x~ v j v+n {x; q), 

fulfil the required conditions. 

We will now do a perturbation analysis of these minimal solutions in a similar way as 
is done for orthogonal polynomials [12]. Again we write the recurrence relation as 

(7.9) h± +ltV - (x' 1 + x)h± v (x; q) + h^_ 1>v (x; q) = -xq" +n h±„(x; q), 

and look at this equation as a non-homogeneous second order recurrence relation with non- 
homogeneous term —xq iy+n h^ 1/ (x; q). The homogeneous equation has two simple solutions, 
h±(x; 0) = x^ n . We solve the non-homogeneous recurrence relation using Green functions, 
but now the Green function Gi{n, m) is the solution of 



J n,m 



(7.10) G 2 (n + 1, m) - (x _1 + x)G 2 {n, to) + G 2 (n - 1, to) = 5 r , 
with boundary conditions 

(7.11) G 2 (n,m) = 0, n < m. 

Since G 2 {m, to) = G 2 {m — 1, to) = we find G 2 {m + 1, to) = 1 and in general 

x + x~ ! 



G 2 (n, to) = Un-m-x y J, n>TO. 

Multiply the recurrence (7.9) by G 2 (n, to) and (7.10) by h^ l/ (x; q) and subtract to find 

h n+iA x ' ?) G '2(n, to) - h* v (x; q)G 2 (n - 1, to) + h^(x; q)5 n , m 

= h* v (x; q)G 2 (n + 1, to) - h^_ l u (x; q)G 2 (n, m) - xq u+n h±^(x; q)G 2 (n, to). 
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Add the equations from n = m to n = M, with m < M and use the boundary conditions 
(7.11) to find 



h %+i,u( x > ^(M, m) - h± tV {x; q)G 2 (M + 1, to) 

M 
n=m+l 

From (7.8) together with (3.7) we obtain 

lim h+j +lu (x;q)G 2 (M,m) - v (x;q)G 2 (M + l,m) = -x" m , |x| > 1, 

and 

lim h- M+lv {x;q)G 2 {M,m) - u (x;q)G 2 (M + l,m) - x m , \x\ < 1, 
so by letting M — > oo we have 

oo / I -1 \ 

(7-12) h± u (x;q) = x^-x £ (s; ? )tf fc _ n _ x ) ■ 

k=n+l ^ ' 

Compare these relations to (7.5). We can find appropriate bounds on these solutions and 
from this we can obtain bounds for J l/+n (x~ 1 ; q) and j l/+n (x; q). 



Lemma 7.2. If x ^ ±1 then 



(7.13) 



and 



(7.14) 



\x n h+ v (x;q)\<exp ^ l _ x _ 2 ^ — - J 

/ 2 g H-n+l\ 



1 _ q + (1-q)2 
!_ g + ( 1 _ Q )2 j 



M > 1, 
kl < 1, 



M > 1, 

Id < 1. 



Proof. We now use a backward version of GronwalPs inequality: for non- negative A, c n , d r 
(n > 0) we have 

oo / oo \ 

c n < A + ^ rfjtCjt < oo c n < Aexp ^ d fe J . 

fe=n+l \fc=n+l / 
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The inequalities (7.13) then follow from (7.12) and the inequalities \x ±n U n ((x+x~ 1 )/2) | < 
2/| 1 — x ±2 \, which hold for |x| < 1 (for the + sign) and \x\ > 1 (for the — sign). 

Inequality (7.14) uses the inequality \x ±n U n {{x + x _1 )/2)| < n + 1 on \x\ < 1 and 
\x\ > 1 respectively. So from (7.12) we get 

oo 

\x n h+ v (x;q)\<l+ q u+k \x k h+ u (x-q)\(k-n) 

k=n+l 

oo 

<1 + kq» +k \x k h+ v (x;q)\, \x\ > 1, 

k=n+l 

from which the first inequality of (7.14) follows by Gronwall's inequality. □ 

We are now ready to give some information about the zeros of the functions h^ u (x; q) 
inside and outside the open unit disk. 

Theorem 7.3. The zeros ofh^ u (x;q) are all real. The function h^ v (x;q) has no zeros 
in {x E C : \x\ > 1} and h~ v (x;q) has no zeros inside {x G C : |x| < 1}, whenever 
n > M(v, q), where 

(7.15) M(v,q) = -v-l + ^ 1 - q) 1 



In q Inq 

In particular hZi u (x; q) has at most 2M(u, q)+2 zeros in {x e C : |x| > 1} and hZi q) 
has at most 2M(v, q) + 2 zeros in {x e C : |x| < 1}. 

Proof. The reality of the zeros follows from the explicit representation (7.8) and the re- 
ality of the zeros of the Hahn-Exton Bessel function [17, §3] and the zeros of j v (x;q), cf. 
theorem 5.3. For an upper bound on the number of zeros, we use (7.12) to find 

1 - q) = £ q» +k x± k h±„(x; q)x^ k ^U M ^f- . 



k=n+l 

Use the inequality (7.14) and \x ±n U n ({x + x~ " 1 )/2) | < n + 1 to find for |x| > 1 

7 i/+fc+l n u+k+l 

k-n)q» +k exp[- 

k=n+l 

and similarly for \x\ < 1 



\i-z n h+ v (x ]q )\< £ ( k - n W +k ^p{-\z^ + fr^w) 



1 — x~ 



l h- v {x;q)\< ^(fc-nrexpU— + 

k=n+i v q V 1) / 
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The right hand side can be bounded by 

°° / U n v+k+l n u+k+l \ / n v+l \ 00 

JJ (*- n)9 ~»«p H— + « < exp « E (,_ nV+ * 

fc=n+l V y \ HJ / \\ H) / k=n+1 

Choose M = M(u,q) such that 

then h^ l/ (x;q) for n > M cannot be zero for any x such that |x| > 1. An appropriate 
M(y,q) is given by (7.15). The same reasoning holds for h~„(x;q) on the closed unit 
disk. So now we have established that for n > M the function fo+ has no zeros for 
\x\ > 1 and h~ v has no zeros for \x\ < 1. The zeros of h+ u are equal to the zeros of 
J u+n (l/x; q). If k = 1, 2, 3, ... , are the zeros of J v {x\ q) numbered increasingly, then 
from the interlacing property of theorem 3.7 in [17] we have j% < < hence when 
the parameter v is decreased by one, then the /cth positive zero moves to the left. This 
means that the /cth positive zero (counted from the right) of h^_ 1 v {x\ q) is to the right of 
the kth positive zero of h^ u (x;q). Since h~l Il/ (x;q) has no zeros x > 1, this means that 
h~M_i v {x]q) can have one zero x > 1, namely l/j] /+M_1 , and it cannot have two zeros 
x > 1 since l/j2 +M_1 < l/j\ +M < 1. Decreasing the degree of h^ l/ (x;q) by one thus 
increases the number of zeros in \x\ > 1 by at most 2 (one positive zero and one negative 
zero). Therefore h't 1 ^(x;q) has at most 2M + 2 zeros in \x\ > 1. A similar reasoning 
works for the zeros of h~ „(x; q) in \x\ < 1 by using the interlacing property of the zeros of 
j u (x; q) and j u+ i(x; q) given by theorem 5.3. □ 

The upper bound on the number of zeros of u (x; q) gives a useful upper bound on the 
number of discrete mass points of the strong moment functional C as given in theorem 6.3 
when s = 1. Indeed, the zeros of h_ lv (x;q) correspond with the zeros of J v -i(l/x; q) 
and thus < M(v, q) + 1. Similarly the zeros of hZi v (x] q) correspond with the zeros of 
j„-i(x; q) and thus M < M{v, q) + 1. In particular, M = N = in theorem 6.3 for s = 1 
for v satisfying M(v, q) < 0. 

Finally let us give another derivation of the orthogonality of the Laurent polynomials 
h n , v {x;q) by using the minimal solutions v (x; q) . Observe that from (3.3), (3.5) and 
(7.8) it follows that h^ v (x; q) have a power expansion of the form 

00 

x n h+ v (x; q) = l + J2 k ) x ~ 2k ' \ x \ > !> 

fc=0 

and 

00 

x~ n h~ jV (x;q) = 1 + ^2K~(n,k)x 2k , \x\ < 1. 
fc=i 
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We can get some information on the coefficients K ± (n, k) by introducing Banach algebras. 
If / is analytic in the open unit disk with Taylor series 

oo 

then we define 

oo 

u/ii-=i>*i/*i, 

k=0 

and we denote by A~ all the functions / for which ||/||_ < oo. Here v^, k G Z + , is a 
positive increasing sequence for which vq = 1 and v n < u m u n - m for every n > m > 0. 
Similarly, when g is analytic near infinity with Laurent series 

oo 
k=0 

then we define 

oo 

= ^2 U k\9k\, 

k=0 

and denote by A + all the functions g for which ||^||+ < oo. One easily verifies that for two 
functions /i , /2 G A ± one has 

||/l/2||±< IIAII± \\h\\±, 

so that we are dealing with Banach algebras. 
Observe that 

oo 

\\x n h+ jl/ (x;q)\\ + = 1 + J2"2k\K+(n,k)l 

fc=0 
oo 

\\x- n h- v (x;q)\\- = l + J2^k\K-(n,k)\. 

k=i 

Taking norms in (7.12) gives 

oo 

\\x n h+ v {x;q)\\ + <l+ q u+k \\x k ht^Q)\\ + \\^- k+1 Uk-n-i 

k=n+l 

Now 

(, -1 \ n n 

X -^) 11+ = II E^" 2J H+ = E^' < (n+ IK 
' 3=0 3=0 
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so that Gronwall's inequality gives 

\\x n hlAx;q)\\+<e W ( jr kq» +k v 2k \. 

\k=n+l J 

Taking v n = a n with a < q~ x l 2 shows that 

oo 

l + ^a 2fc |K+(n,fc)| < oo, 

k=0 

so that x n h+ v {x; q) G A + . This shows that the function h+ v {x;q) is in fact defined for 
x > q 1 ! 2 . A similar reasoning shows that x~ n h~ J/ (x; q) £ A~ and that h~ v (x; q) is defined 
for x < q~ 1 / 2 . From (7.8) we see that h+ u (x;q) has poles at the zeros of (qx~ 2 ; q)^ and 
that x = ±q x l 2 are the poles of largest modulus. Similarly h~ v {x; q) has poles at the zeros 
of {qx 2 ; 5)00 and x = ±(/ -1 / 2 are the poles of smallest modulus. 

Suppose now that ±1 are not zeros of h^ 1 (x; q). Evaluate the contour integral 

(x; q)h+ u (x; q) . 

= / t — - — L - ax. 

2mJ\ x \ =1 hZx Jx;q) 



If m < n then near x = 00 the integrand behaves as x m n 1 and thus I + has no contri- 
bution from x = 00. So when x^ (j > 1) are the zeros of hZi v (x; q), then 



N 



h m ,v(xj;q)h+ (xj;q) 



[ht ltV (x+;q)]' ' 

where N is defined as in theorem 6.3 for s = 1. Similarly we compute the contour integral 



1 f h m ^ l/ (x;q)h n>u (x;q) 
2ni J\x\=i hZ x Ax]q) 



The integrand behaves as x n ~ m+1 near x = and thus there is no pole at the origin when 
m < n (even for m < n + 1). There are poles at the zeros xj (j > 1) of hZ± >u (x; q) and 
we thus have 

I _ h m , v {xj;q)h-^{xj;q) 
~~j±L [hZ ltV (xJ ]q )]' ' 

where M is defined as in theorem 6.3 for s = 1. Subtracting I + and 7_ gives 

r T 1 /" , , v fKA x '> < fi h -i,»( x > < fi ~ h nA X '^) h -lA X '^) 1 . 

J - - J + = 77- / ^m,^(^; g) r+ — ; — rr= — — ; dx. 

2m J\x\=i V 
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The Laurent polynomial h njU (x; q) is a solution of the three-term recurrence relation (7.1) 
and therefore a linear combination of the two special solutions h^ u (x; q). With the initial 
conditions ho jU (x;q) = 1 and h-i- u (x; q) = and by combining (7.8) with lemma 6.2 we 
find 

Or" 1 - x)hn, v (x; q) = hr lt „(x; q)h+ v (x; q) - h\„(x; q)h~ v (x; q), 

so that 

I--I+ = ^— [ (x - x- 1 )h n ^(x;q)h m!U (x;q)- T — — -. 

lm J\x\=i hZ hl/ (x; q)h_ l v {x] q) 

On the other hand, at a zero x^ we see that hZjy(Xj~;q) is a solution of (7.1) with initial 
value h_ l u (Xj~;q) = 0, so that hZ jU (x^;q) = fiQ l ,(Xj,q)h n!l/ (Xj,q). Similarly at a zero 
xj we have h~ v (x~\q) = h^ v {x~ ; q)h n>v (xj ; q). Therefore 

M{v,q) h + u (x~;q) 



[^i>7;?)]' 

M{v,q) h + (x + -q) 

Combining both expressions for 7_ — / + gives the orthogonality of the Laurent polynomials 
h n ,iy(x] q) and corresponds to the result given in theorem 6.3 for s = 1. In case q 1 / 2 < s < 
g -1 / 2 the orthogonality relations of theorem 6.3 can be derived in a similar way. 

This approach can also be used to prove the orthogonality for the Laurent polynomials 
x~ x h n)U {x; q), cf. (2.3). Note also that the case q = gives the orthogonality relations for 
the Chebyshev polynomials of the second kind. 
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